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Abstract 

We investigate the environments of quasars such as number distribution of galaxies 
using a semi-analytic model which includes both galaxy and quasar formations based 
on the hierarchical clustering scenario. We assume that a supermassive black hole 
is fueled by accretion of cold gas and that it is a source of quasar activity during 
a major merger of the quasar host galaxy with another galaxy. This major merger 
causes spheroid formation of the host galaxy. Our model can reproduce not only 
general form of the galaxy luminosity functions in the local Universe but also the 
observed relation between a supermassive black hole mass and a spheroid luminosity, 
the present black hole mass function and the quasar luminosity functions at different 
redshifts. Using this model, we predict the mean number of quasars per halo, bias 
parameter of quasars and the probability distribution of the number of galaxies around 
quasars. In our model, analysis of the mean number of quasars per halo shows that 
the spatial distribution of galaxies is different from that of quasars. Furthermore, 
we found from calculation of the probability distribution of galaxy numbers that at 
0.2 < 2 < 0.5, most quasars are likely to reside in galaxy groups. On the other hand, 

1 ^ ^ ^ 2 most quasars seem to reside in more varied environments than at a 
lower redshift; quasars reside in environments ranging from small groups of galaxies 
to clusters of galaxies. Comparing these predictions with observations in future will 
enable us to constrain our quasar formation model. 

Key words: galaxies: evolution — galaxies: formation — galaxies: quasars: 
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1. Introduction 



The environments of quasars provide important clues to the physical processes of their 
formation and also yield important information about the relations between the distribution 
of quasars and the large-scale structure of the universe. For more than three decades, we have 
known that quasars are associated with enhancements in the spatial distributions of galaxies 
(Bahcall, Schmidt, Gunn 1969). Studies of the environments of quasars in the nearby universe 
(z ^ 0.4) have shown that quasars reside in environments ranging from small to moderate groups 
of galaxies rather than in rich clusters (e.g. Bahcall, Chokshi 1991; Fisher et al. 1996; McLure, 
Dunlop 2001). In order to interpret the observational results of the environments of quasars 
at low redshifts and predict the environments of quasars at high redshifts, a physical model of 
quasar formation based on cosmological context is required. 

It has become widely accepted that quasars are fueled by accretion of gas onto super- 
massive black holes (SMBHs) in the nuclei of host galaxies since Lynden-Bell (1969) proposed 
this idea on quasars. Recent observations of galactic centers suggest that a lot of nearby 
galaxies have central black holes and their estimated masses correlate with the luminosities 
of spheroids^ of their host galaxies (e.g. Kormendy, Richstone 1995; Magorrian et al. 1998; 
Merritt, Ferrarese 2001). The connection between SMBHs and their host spheroids suggests 
that the formation of SMBHs physically links the formation of the spheroids which harbor the 
SMBHs. Thus, this implies that the formation of quasars is closely related to the formation of 
galaxies, especially of spheroids. Therefore, in order to study the formation and evolution of 
quasars, it is necessary to construct a unified model which includes both galaxy formation and 
quasar formation. 

Recently, some authors have tried to construct galaxy formation models on the basis of 
the theory of hierarchical structure formation in cold dark matter (CDM) universe. These efforts 
are referred to as semi-analytic models (SAMs) of galaxy formation. In the CDM universe, 
dark matter halos cluster gravitationally and merge together in a manner that depends on 
the adopted power spectrum of initial density fluctuations. In each of the merged dark halos, 
radiative gas cooling, star formation, and supernova feedback occur. The cooled dense gas 
and stars constitute galaxies. These galaxies sometimes merge together in a common dark 
halo and more massive galaxies form. In SAMs, the merger trees of dark matter halos are 
constructed using a Monte-Carlo algorithm and simple models are adopted to describe the 
above gas processes. Stellar population synthesis models are used to calculate the luminosities 
and colors of model galaxies. It is therefore straightforward to understand how galaxies form 
and evolve within the context of this model. SAMs successfully have reproduced a variety of 
observed features of local galaxies such as their luminosity functions, color distribution, and so 
on (e.g. Kauffmann, White, Guiderdoni 1993; Cole et al. 1994, 2000; Somerville, Primack 1999; 

^ Throughout this paper, we refer to bulge or eUiptical galaxy as spheroid. 
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Nagashima et al. 2001, 2002). 

In these models, it is assumed that disk stars are formed by coohng of gas in the halo. 
If two galaxies of comparable mass merge, it is assumed that starbursts occur and form the 
spheroidal component in the center of the galaxy. A^-body simulations have shown that a 
merger hypothesis for the origin of spheroids can explain their detailed internal structure (e.g. 
Barnes 1988; Hernquist 1992, 1993; Heyl, Hernquist, Spergel 1994). Kauffmann and Chariot 
(1998) have demonstrated that the merger scenario for the formation of elliptical galaxies is 
consistent with the color-magnitude relation and its redshift evolution (see also Nagashima, 
Gouda 2001). On the other hand, hydrodynamical simulations have shown that a merger 
of galaxies drives gas to fall rapidly to the center of a merged system and to fuel nuclear 
starburst (Negroponte, White 1983; Mihos, Hernquist 1994, 1996; Barnes, Hernquist 1996). 
Moreover, observed images of quasar hosts show that many quasars reside in interacting systems 
or elliptical galaxies (Bahcall et al. 1997). Therefore, it has often been thought that the major 
merger of galaxies would be a possible mechanism for quasar and spheroid formation. 

So far, a lot of studies on quasar evolution based on the hierarchical clustering scenario 
have been carried out with the assumption that the formation of quasars is linked to the first 
collapse of dark matter halos with galactic mass and that these models can explain the decline 
of quasar number density at 2; ^ 3 (e.g. Efstathiou, Rees 1988; Haehnelt, Rees 1993) and 
properties of luminosity functions of quasars (e.g. Haiman, Loeb 1998; Haehnelt, Natarajan, 
Rees 1998; Hosokawa et al. 2001). However, if quasars are directly linked to spheroids of host 
galaxies rather than to dark matter halos, the approximation of a one-to-one relation between 
quasar hosts and dark matter halos would be very crude, especially at low redshift. Therefore, 
it is necessary to construct a model related to spheroid formation and SMBH formation directly. 
Kauffmann and Haehnelt (2000) introduced a unified model of the evolution of galaxies and 
quasars within the framework of SAM (see also Cattaneo 2001). They assumed that SMBHs 
are formed and fueled during major galaxy mergers and their model reproduces quantitatively 
the observed relation between spheroid luminosity and black hole mass in nearby galaxies, 
the strong evolution of the quasar population with redshift, and the relation between the 
luminosities of nearby quasars and those of their host galaxies. 

In this paper, we investigate properties of quasar environments, using a SAM incorpo- 
rated simple quasar evolution model. We assume that SMBHs are formed and fueled during 
major galaxy mergers and the fueling process leads quasar activity. While this assumption 
is similar to the model of Kauffmann and Haehnelt (2000), our galaxy formation model and 
the adopted model of fueling process are different from their model. Here we focus on optical 
properties of quasars and attempt to consider the number of quasars per halo, effective bias 
parameter of quasars and the number of galaxies around quasars as characterizations of en- 
vironments of quasars, because a) these quantities provide a direct measure of bias in their 
distribution with respect to galaxies and b) comparing results of the model with observations 
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will enable us to constrain our quasar formation model. 

The paper is organized as follows: in §2 we briefly review our SAM for galaxy forma- 
tion; in §3 we introduce the quasar formation model; in §4 we calculate the galaxy number 
distribution function around quasars; in §5 we provide a summary and discussion. 

In this study, we use a low-density, spatially fiat cold dark matter (ACDM) universe 
with the present density parameter f2o = 0.3, the cosmological constant Aq = 0.7, the Hubble 
constant in units of 100km s~^ Mpc~^ h = 0.7 and the present rms density fluctuation in spheres 
of 8/i~^Mpc radius as = 1.0. 

2. Model of Galaxy Formation 

In this section we briefly describe our SAM for the galaxy formation model, details of 
which are shown in Nagashima et al. (2001). Our present SAM analysis obtains essentially the 
same results as the previous SAM analyses, with minor differences in a number of details. 

First, we construct Monte Carlo realizations of merging histories of dark matter halos 
using the method of Somerville, Kolatt (1999), which is based on the extended Press-Schechter 
formalism (Press, Schechter 1974; Bower 1991; Bond et al. 1991; Lacey, Cole 1993). We adopt 
the power spectrum for the specific cosmological model from Bardeen et al. (1986). Halos with 
circular velocity Vchc <40km s~^ are treated as diffuse accretion matter. The evolution of the 
baryonic component is followed until the output redshift coincides with the redshift interval 
of Az = 0.06(1 + z), corresponding to the dynamical time scale of halos which collapse at the 
redshift z. Note that Shimizu et al. (2002) recently pointed out that a much shorter timestep 
is required to correctly reproduce the mass function given by the Press-Schechter formalism. 
However, a serious problem exists only at small mass scales (^ lO^^M©). Thus we use the above 
prescription of timestep. 

If a dark matter halo has no progenitor halos, the mass fraction of the gas in the halo 
is given by Q-^/Qq , where Qh = 0.015/i~^ is the baryonic density parameter constrained by 
primordial nucleosynthesis calculations (e.g. Suzuki, Yoshii, Beers 2000). Note that a recent 
measurement of the anisotropy of the cosmic microwave background by the BOOMERANG 
project suggests a slight higher value, l^b ~ 0.02/i"^ (Netterfield et al. 2002). Cole et al. (2000) 
have already investigated the effect of changing f2b and showed that this mainly affects the 
value of the invisible stellar mass fraction such as brown dwarfs parameterized by T (see 
below). When a dark matter halo collapses, the gas in the halo is shock-heated to the virial 
temperature of the halo. We refer to this heated gas as the hot gas. At the same time, the 
gas in dense regions of the halo cools due to efficient radiative cooling. We call this cooled gas 
the cold gas. Assuming a singular isothermal density distribution of the hot gas and using the 
metallicity-dependent cooling function by Sutherland, Dopita (1993), we calculate the amount 
of cold gas which eventually falls onto a central galaxy in the halo. In order to avoid the 
formation of unphysically large galaxies, the above cooling process is applied only to halos with 
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Vcirc <400 km s~^. This handling would be needed because the simple isothermal distribution 
forms so-called "monster galaxies" due to the too efficient cooling at the center of halos. While 
this problem will probably solved by adopting another isothermal distribution with central core 
(Cole et al. 2000), we take the above approach for simplicity. 

Stars are formed from the cold gas at a rate of = Mcom/t*, where Mcom is the mass 
of cold gas and r* is the time scale of star formation. We assume that is independent of z, 
but dependent on V^irc as follows: 

The free parameters of r° and are fixed by matching the observed mass fraction of cold gas 
in neutral form in the disks of spiral galaxies. In our SAM, stars with masses larger than lOM© 
explode as Type II supernovae (SNe) and heat up the surrounding cold gas. This SN feedback 
reheats the cold gas to hot gas at a rate of Mreheat = /3M*, where (3 is the efficiency of reheating. 
We assume that (3 depends on Vcirc as follows: 
/y. \-"hot 

P={^) ■ (2) 

The free parameters of Vhot and tthot are determined by matching the local luminosity function 
of galaxies. With these and Mreheat thus determined, we obtain the masses of hot gas, cold 
gas, and disk stars as a function of time during the evolution of galaxies. 

Given the star formation rate as a function of time, the absolute luminosity and colors 
of individual galaxies are calculated using a population synthesis code by Kodama, Arimoto 
(1997). The initial stellar mass function (IMF) that we adopt is the power- law IMF of Salpeter 
form with lower and upper mass limits of O.IM0 and GOMq, respectively. Since our knowledge 
of the lower mass limit is incomplete, there is the possibility that many brown dwarf-like 
objects are formed. Therefore, following Cole et al. (1994), we introduce a parameter defined 
as T = (Mium + MBD)/Mium, where Mium is the total mass of luminous stars with m > O.IMq 
and Mbd is that of invisible brown dwarfs. To account for extinction by internal dust we adopt 
a simple model by Wang, Heckman (1996) in which the optical depth in i?-band is related to 
the luminosity as tb = 0.8(^^/1.3 x lO^^L©)"'^. Optical depths in other bands are calculated 
by using the Galactic extinction curve, and the dust distribution in disks is assumed to be the 
slab model considered by Somerville, Primack (1999). 

When several progenitor halos have merged, the newly formed larger halo should contain 
at least two or more galaxies which had originally resided in the individual progenitor halos. 
We identify the central galaxy in the new common halo with the central galaxy contained in the 
most massive of the progenitor halos. Other galaxies are regarded as satellite galaxies. These 
satellites merge by either dynamical friction or random collision. The time scale of merging by 
dynamical friction is given by 



5 




""'"^ " InA, Up J ll03kms-vil0i2Moj 
where i?H and Kirc are the radius and the circular velocity of the new common halo, respectively, 
In Ac is the Coulomb logarithm, and Mgat is the mass of the satellite galaxy including its dark 
matter halo (Binney, Tremaine 1987). When the time passed after a galaxy becomes a satellite 
exceeds Tfj-ic, a satellite galaxy infalls onto the central galaxy. On the other hand, the mean free 
time scale of random collision is given by 

^2 l^Mpc ; 1^0.12Mpc J 

X ( "^-^ , ) ~" ( " Gyr, (4) 

VlOOkms^V VaOOkms^V ^' 

where is the number of satellite galaxies, rgai is their radius, and dhaio and (jgai are the ID 
velocity dispersions of the common halo and satellite galaxies, respectively (Makino, Hut 1997). 
With a probability of At/rcou, where At is the timestep corresponding to the redshift interval 
A2;, a satellite galaxy merges with another randomly picked satellite. 

Consider the case that two galaxies of masses mi and m2(> mi) merge together. If the 
mass ratio / = mi/m2 is larger than a certain critical value of /bulge, we assume that a starburst 
occurs and all the cold gas turns into stars and hot gas, which fills the halo, and all of the stars 
populate the bulge of a new galaxy. On the other hand, if / < /bulge, no starburst occurs and a 
smaller galaxy is simply absorbed into the disk of a larger galaxy. These processes are repeated 
until the output redshift. We classify galaxies into different morphological types according 
to the i?-band bulge-to-disk luminosity ratio B/D. In this paper, galaxies with B/D > 2/3, 
and B/D < 2/3 are classified as ellipticals/SOs and spirals, respectively. This method of type 
classification well reproduces the observed type mix. 

The above procedure is a standard one in the SAM for galaxy formation. Model param- 
eters are determined by comparison with observations of the local Universe. In this study, we 
use the astrophysical parameters determined by Nagashima et al. (2001) from local observa- 
tions such as luminosity functions, and galaxy number counts in the Hubble Deep Field. The 
adopted parameters of this model are tabulated in Table 1. In Figure 1 we plot the results 
of local luminosity functions of galaxies represented by solid lines. Note that the resultant 
luminosity functions hardly change if the SMBH formation model is included (dashed lines; 
see the next section). Symbols with errorbars indicate observational results from the 5-band 
redshift surveys (APM, Loveday et al. 1992; 2dF, Folkes et al. 1999) and from the /C-band 
redshift surveys (Gardner et al. 1997; 2MASS, Cole et al. 2001). As can be seen, the results of 
our model using these parameters are generally consistent with the observations, at least with 
the APM result. 
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Table 1. Model Parameters 



cosmological parameters 




astrophysical parameters 


flo Ao h fTg 


Vhot (km s" 


"hot T° (Gyr) /bulge T 


0.3 0.7 0.7 1 


280 


2.5 1.5 -2 0.5 1.5 




Mb-5 log[h] Mk-5 log[h] 

Fig. 1. Local luminosity fimctions in the (a) B-band and (b) iiT-band. The thick line shows the result 
of the model without SMBH formation. The short dashed line shows the model with SMBH formation. 
Symbols with errorbars in (a) indicate the observational data from APM (Loveday et al. 1992, filled 
squares) and 2dF (Folkes et al. 1999, open circles). Symbols in (b) indicate the data from Gardner et al. 
(1997; open inverted triangles), and 2MASS (Cole et al. 2001, filled triangles). 

3. Model of Quasar Formation 

In this section, we introduce a quasar formation and evolution model into our SAM. As 
mentioned earlier, the masses of SMBHs have tight correlation with the spheroid masses of their 
host galaxies (e.g. Kormendy, Richstone 1995; Magorrian et al. 1998; Merritt, Ferrarese 2001) 
and the hosts of quasars found in the local Universe are giant elliptical galaxies or galaxies 
displaying evidence of major mergers of galaxies (Bahcall et al. 1997). Moreover, in SAMs for 
galaxy formation, it is assumed that a galaxy-galaxy major merger leads to the formation of 
a spheroid. Therefore, we assume that SMBHs grow by merging and are fueled by accreted 
cold gas during major mergers of galaxies. When host galaxies merge, pre-existing SMBHs sink 
to the center of the new merged galaxy owing to dynamical friction and finally coalesce. The 
timescale for this process is unknown, but for the sake of simplicity we assume that SMBHs 
merge instantaneously. Gas-dynamical simulations have demonstrated that the major merger 
of galaxies can drive substantial gaseous infiows and trigger starburst activity (Negroponte, 
White 1983; Mihos, Hernquist 1994, 1996; Barnes, Hernquist 1996). Thus, we assume that 
during major merger, some fraction of the cold gas that is proportional to the total mass of stars 
newly formed at starburst is accreted onto the newly formed SMBH. Under this assumption, 
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the mass of cold gas accreted on a SMBH is given by 



acc 



/bhAM, 

,burst ) 



(5) 



where /bh is a constant and AM*_burst is the total mass of stars formed at starburst. AM*^burst 
is derived in the Appendix. The free parameter of /bh is fixed by matching the observed 
relation between a spheroid luminosity and a black hole mass found by Magorrian et al. (1998) 
and we find that the favorable value of /bh is nearly 0.03. In Figure 2 we show scatterplots 
(open circles) of the absolute ^-band magnitudes of spheroids versus masses of SMBHs of 
model for /bh = 0.03. The thick solid line is the observational relation and the dashed lines 
are the la scatter in the observations obtained by Magorrian et al. (1998). For /bh ^ 0.3, 
changing /bh shifts the black hole mass almost linearly. The obtained gas fraction (/bh = 0.03) 
is so small that the inclusion of SMBH formation does not change the properties of galaxies 
in the local Universe. In Figure 1, the dashed lines show the results of the model with the 
SMBH formation. This result differs negligibly from the result of the model without SMBH 
formation. Therefore, we use the same astrophysical parameters tabulated in Table 1 regardless 
of inclusion of the SMBH formation model. Figure 3 (a) shows black hole mass functions in our 
model at a series of redshifts. This indicates that the number density of the most massive black 
holes increases monotonically with time in the scenario where SMBHs grow by accretion of gas 
and by merging. In Figure 3 (a), we superpose the present black hole mass function obtained 
by Salucci et al. (1999). They derived this black hole mass function from the observed radio 
luminosity function of nearby radio-quiet galaxies and the empirical correlation between radio 
luminosities from the nuclei of radio-quiet galaxies and the mass of their black holes. Our model 
result is consistent with their mass function. For comparison, we also plot the mass functions 
of bulge and disk for all galaxies in Figure 3 (b) and (c), respectively. The steep slopes at low 
masses of mass functions of black hole and bulge are mainly due to random collisions between 
satellite galaxies in this model. 

To obtain the observed linear relation between a spheroid luminosity and a black hole 
mass, Kauffmann and Haehnelt (2000) adopted model of fueling process in which the ratio of 
accreted mass to total available cold gas mass scales with halo circular velocity in the same way 
as the mass of stars formed per unit mass of cooling gas. While their approach is similar to 
ours, their star formation and feedback models are different from ours and they do not consider 
random collisions. Therefore, their resultant model description is slightly differ from ours in 
equation (5). 

Next, we consider the light curve of quasars. We assume that a fixed fraction of the rest 
mass energy of the accreted gas is radiated in the i?-band and the quasar life timescale tiife(^) 
scales with the dynamical time scale tdyn of the host galaxy where tdyn oc rgai/cigai oc Rn/Vdrc- 
Here we adopt the 5-band luminosity of a quasar at time t after the major merger as follows; 



Lsit) = LB(peak)exp(-t/tiifc). 



(6) 
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Fig. 2. The relation absolute y-band spheroid magnitude - mass of SMBH. The open circles are an abso- 
lute y-band magnitude limited sample of spheroids in our model. The thick solid line is the observational 
relation obtained by Magorrian et al. (1998). The dashed lines indicate the Icr scatter in the observations. 

The peak luminosity -^^(peak) is given by 

Lb (peak) = ££^^^££^^ (7) 

^life 

where is the radiative efficiency in S-band, tufe is the quasar hfe timescale and c is the speed 
of hght. In order to determine the parameter €3 and the present quasar hfe timescale tiifc(O), 
we have chosen them to match our model luminosity function with the observed abundance 
of bright quasars at z = 2. We obtain = 0.005 and tiifc(O) = 3.0 x lO^yr. The resulting 
luminosity functions at four different redshifts are shown in Figure 4. We superpose the lumi- 
nosity functions derived from the 2dF 10k catalogue (Croom et al. 2001) for a cosmology with 
Qq = 0.3, Ao = 0.7 and h = 0.7, which is analyzed and kindly provided by T. T. Takeuchi. He 
used the method of Efstathiou, Ellis, Peterson (1988) for the estimation of the luminosity func- 
tions. In order to reanalyze the error with greater accuracy, they applied bootstrap resampling 
according to the method of Takeuchi, Yoshikawa, Ishii (2000). Absolute S-band magnitudes 
were derived for the quasars using the fc-corrections derived by Cristiani, Vio (1990). Our model 
reproduces reasonably well the evolution of observed luminosity functions. Thus, in the next 
section, we use these model parameters in order to investigate the environments of quasars. 

For comparison, we also plot the result of model with = 0.005 and tiifc(O) = 3.0 x lO^yr 
in Figure 4 (dot-dashed lines). In this case, the abundance of luminous quasars decreases. To 
prolong a quasar life timescale affects the quasar luminosity function due to the following two 
factors: a decrease in the peak luminosity Lb (eq.[7]) and an increase in the exponential factor 
exp(— t/tiife) in equation (6). For the majority of bright quasars, the elapsed time t since the 
major merger is much smaller than the quasar life timescale tufc, t/tufc <^ 1- Therefore, the 
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Fig. 3. (a) The black holes mass function of models for /bh = 0.03 as a function of epoch. The solid, 
short-dashed and dot-dashed lines indicate the results at z = 0, 1 and 2 respectively. The symbols with 
errorbars are the present black hole mass function obtained by Salucci et al. (1999). (b) The bulge and (c) 
the disk mass functions of model galaxies as a function of epoch. The solid, short-dashed and dot-dashed 
lines indicate the results at z = 0, 1 and 2 respectively. 
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Fig. 4. The i?-band quasar luminosity functions at (a)z = 0.25, (b) z = 0.5, (c)z = 1.25 and (d)z = 2.0. The 
sohd hues are iiifc(O) = 3 x lO'^yr and the dot-dashed hnes are tiifo(O) = 3 x lO^yr. The symbols show results 
from the 2dF 10k catalogue (Croom et al. 2001) reanalyzed by Takeuchi for a cosmology JIq — 0.3, Xq — 0.7 
and ft, = 0.7. 

former factor dominates the latter and the number of luminous quasars decreases. Thus, a 
long quasar life timescale results in a very steep quasar luminosity function. Note that if we 
change the radiative efficiency e^, the quasar luminosities simply scale by a constant factor in 
our model. Thus, changing shifts the luminosity function horizontally. 

4. Environments of Quasars 

In this section, we investigate the environments of quasars using our model. We con- 
sider the halo mass dependence of the mean number of quasars per halo and the probability 
distribution of the number of galaxies around quasars as characterizations of the environments 
of quasars. This is because the former is one of measures of the relation between quasars and 
dark matter distributions and the latter reflects the relationship between galaxies and quasars. 

In Figure 5, we plot (Ngsi^M)) and {Nqso{M)) that denote the mean number of galaxies 
and quasars per halo with mass M, respectively, at (a) z = 0.5 and (b) z = 2.0. We select galaxies 
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with Mb — 5\og{h) < — 19 and quasars with Mb — 51og(/i) < —21, where Mb is absolute 5-band 
magnitude. It should be noted that changing the magnitude of selection criteria for galaxies 
and quasars would alter these results, but qualitative features are not altered. As is seen in 
Figure 5, there are more galaxies and quasars at high z. At higher redshift, halos have more 
cold gas available to form stars and to fuel SMBHs because there has been relatively little time 
for star formation to deplete the cold gas at these redshifts. Thus, the number of luminous 
galaxies grows. Furthermore, at higher redshift, both timescales of the dynamical friction and 
the random collisions are shorter because the mass density of a halo is higher. Therefore, the 
galaxy merging rate increases. Consequently, the number of quasars also grows. Moreover, the 
decrease in the quasar life timescale tufe with redshift also contributes to the increase in the 
number of quasars because quasars become brighter as a result of decrease in tufe (eq. [7]). 

From Figure 5, we find that the dependence of {Nqso{M)) on halo mass M is different 
from the dependence of {Nga^^M)). Furthermore, Figure 6 shows that the ratio of {Nqso{M)) 
to (iVgai(M)) varies with redshift and halo mass. Benson et al. (2000) used a combination of 
cosmological A^-body simulation and semi-analytic modeling of galaxy formation and showed 
that the galaxy spatial distribution is sensitive to the efficiency with which galaxies form in 
halos with different mass. Seljak (2000) also obtained the same conclusion using an analytic 
model of galaxy clustering. These results are applicable to the quasar spatial distribution. 
Therefore, our result indicates that the clustering properties of galaxies are not the same as 
those of quasars and that the bias in the spatial distribution of galaxies relative to that of dark 
matter is not the same as the bias in the spatial distribution of quasars. Assumed that biases 
are independent of scale, we can calculate effective biases using the method of Baugh et al 
(1999) as follows; 

/ b{M, z) {N{M, z))n{M\z)dM 



(8) 



J{N{M,z))n{M\z)dM 

where b{M,z) is the bias parameter for dark matter halos of mass M at z, {N{M,z)) denotes the 
mean number of objects (galaxies or quasars) in a halo of mass M at 2; that satisfy the selection 
criteria and n{M\z) is the dark halo mass function at z. Our SAM adopts the Press-Schechter 
mass function which is given by 



n{M\z)dM 




7rMa^{M) 



da{M) 



dM 



exp 



1 m 

2a2(M) 



dM, 



(9) 



where po is the present mean density of the universe, cr(M) is the rms linear density fluctuation 
on the scale M at 2; = and Sc{z) = 6c/D{z). D{z) is the linear growth factor, normalized 
to unity at the present day and 6c is the linear critical density contrast at the collapse epoch. 
Here, we use an approximate formula of 6c for spatially fiat cosmological model (Nakamura & 
Suto 1997). The bias parameter for dark matter halos is given by Jing (1998); 

(0.06-0. 02ncfr) 



b{M,z) 









k'(M) 1 
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Fig. 5. The mean numbers of galaxies witii Mb — 5 log(/i) < —19 (solid lines) and quasars with 
Mb — 51og(/i) < —21 (dot-dashed lines) per halo with mass AI at (a) z — 0.5 and (b) z = 2.0. The 
horizontal dotted Hue marks (iVgai) = 1 and {Nqsq) = 1 

where riefr is the effective spectral index of the power spectrum, din P{k)/ dink, at the wavenum- 
ber defined by the Lagrangian radius of the dark matter halo, k = 2n/rL and r^, = {3M/Aitpq)^^^ . 
Figure 7 shows the evolution of effective bias for galaxies with Mb — 5log{h) < —19 and quasars 
with Mb — 51og(/i) < —21. As is seen in Figure 7, quasars are higher biased tracer than galaxies. 
Furthermore, the evolution of quasar bias is different from that of galaxy bias. This reflects 
the difference in th dependence on halo mass M and redshift of {Nqso{M,z)) and {Nga.i{M,z)) . 
Note that these effective biases are valid for large scale where objects (galaxies or quasars) 
which contribute two-point correlation function populate different halos. 

Next, we formulate the conditional probability that a halo with A'^qso quasars has A'gai 
galaxies. The number density of the halos which contains A'gai galaxies and A^'qso quasars at z 
is obtained from the following expression: 

ri(iVgai,iVQSok) = / N{N,,i,NQso\M,z)n{M\z)dM, (11) 

where A^(A'^gai, A^qsoI^? ^)'^A'gai(iA'^QSO denotes the number of the halos with mass M which 
contains N^^^i ~ A'gai + dN^^^i galaxies and A^'qso ~ -^V'qso + dNqso quasars at z and n{M\z) is the 
dark halo mass function at z. The number density of the halos which contain A'^qso quasars at 
z is obtained from the following expression: 

n(iVQSok) = / N{Nciso\M,z)n{M\z)dM, (12) 

where A^(Ai'Qso |^,^)'^-^qso denotes the number of the halos with mass M which contain A'^qso ~ 
-^Qso + o^-^QSO quasars at z. From equation (11) and (12), the conditional probability that the 
halo with A^'qso quasars has A'gai ~ A^'gai + dNg^^i galaxies at z is given by 

P(iV,a.|iVQSO, ^)rfiV,al = ^^^^f^^ ^^.al ■ (13) 
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log[M/Msij^] 

Fig. 6. The ratio of the mean number of galaxies with Mb — 51og(/i) < — 19 to the mean number of quasars 
with Mb — 51og(/i) < —21 per halo with mass M at z = 0.5 (solid line) and z = 2.0 (dot-dashed line). 




log(l+z) 

Fig. 7. The effective bias parameter of galaxies with Mb — 5\og(h) < — 19 (squares with solid line) and 
quasars with Mb — 51og(/i) < —21 (triangles with dot-dashed line) at z = 0.25, z = 0.5, z = 1.25 and z ^ 2.0. 
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As is seen in the above formulation, given N{Ngsi, Nqso\M, z) and N{Nqso\M, z) from the 
quasar formation model, one can calculate the probability distribution for the number of galax- 
ies around quasars. Figure 8 shows these galaxy number distribution functions around quasars 
estimated by our model. The results are shown for quasars brighter than Mb — 5log{h) = —22 
and for galaxies brighter than Mb — 51og(/i) = —19. Note that at z = 0.25 and z = 0.5 
^'(iVgailiVQSO = 2) = and P{Ng^i\NQso = 3) = for all A^g^i (Fig. 8(a) and (b)) and that 
aX z = 1.25 P{N^iCl\Nqiso = 3) = for all A'gai (Fig. 8(c)). At lower redshift, a halo has at 
most one quasar. Fig 8(a) and (b) show that the halo which has one quasar contains several 
galaxies by high probability. These results indicate that most quasars tend to reside in groups 
of galaxies at 0.2 ^ z ^ 0.5 and is consistent with the observation at z ^ 0.4 (e.g. Bahcall, 
Chokshi 1991; Fisher et al. 1996; McLure, Dunlop 2001). On the other hand, at higher red- 
shift, the numbers of galaxies in the halo with one or two quasars is from several to dozens 
(Fig 8(c) and (d)). These results indicate that quasars locate in ranging from small groups 
of galaxies to clusters of galaxies. Thus at 1 ^ z ^ 2 quasars seem to reside in more varied 
environments than at lower redshift. Kauffmann and Haehnelt (2002) used a combination of 
cosmological iV-body simulation and semi-analytic modeling of galaxy and quasar formation, 
and showed that the ratio of the amplitude of the quasar-galaxy cross correlation function to 
that of the galaxy autocorrelation function decrease with redshift. This indicates that the dif- 
ference between galaxy and quasar distribution becomes smaller at higher redshift. Thus, our 
results obtained by -P(A'gai|Ai'Qso, ^) is not in conflict with their results. 

5. Summary and Discussion 

We have constructed a unified semi-analytic model for galaxy and quasar formation and 
have predicted the mean number of quasars per halo with mass M, {Nqso{M)) , the effective 
bias parameter of quasars &cfr,Qso(^) and probability distribution of the number of galaxies 
around quasars, -P(A'gai|A''Qso), as characterizations of the environments of quasars. These 
quantities reflect the processes of quasar formation such as the amount of cold gas available for 
fueling, the galaxy merger rate and the quasar life timescale. Therefore, by comparing these 
predictions with observations, one will be able to constrain quasar formation models. 

Our model can reproduce not only general form of the galaxy luminosity functions in 
the local Universe but also the observed relation of the SMBH mass to spheroid luminosity, 
and the quasar luminosity functions at different redshifts (Fig.2 and Fig.4). Using this model, 
we have shown (iVQso(M)) and P(A^gai|A^Qso)- The ratio of (iVQso(M)) to (A^gai(M)) varies 
with halo mass in our model (Fig5). These results of our model suggest that the clustering of 
galaxies is not the same as the clustering of quasars and the effective bias parameter of quasars 
and its evolution are different from these of galaxies (Fig. 7). Furthermore, we predict the 
galaxy number distribution function around quasars, -P(Ai'gai|A''Qso) (Fig8). At lower redshifts 
(0.2 ^z^ 0.5), most halos which have quasars have at most several galaxies. This indicates that 
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Fig. 8. The probability distribution for the numbers of galaxies around quasars (a)z = 0.25, (b) z ~ 0.5, 
(0)2: = 1.25 and (d)z = 2.0. The selected galaxies are brighter than Mb — 51og(/i) = —19 and the selected 
quasars are brighter than A^b — 51og(/i) = —22. Long-dashed, solid, dot-dashed and short-dashed lines 
show results for A^qso = 0,1,2 and 3 respectively. 

most quasars reside in groups of galaxies. On the other hand, at higher redshift (1 ^ 2; ^ 2), the 
number of galaxies in the halo with quasars is from several to dozens; quasars reside in ranging 
from small groups of galaxies to clusters of galaxies. These results show that most quasars at 
higher redshift reside in more varied environments than at lower redshift. This model prediction 
is checkable by statistics of galaxies around quasars which will be obtained in future. 

It is still controversial whether the environments of quasars depend on their optical 
and radio luminosities. Some authors have claimed that radio-loud quasars were located in 
richer environments than radio-quiet quasars at at z < 0.6 (e.g. Yee, Green 1984; Yee, Green 
1987; EUingson, Yee, Green 1991; Hintzen, Romanishin, Vlades 1991). However, other people 
obtained a different result. For example, Hutchings, Crampton, Johnson (1995) observed the 
galaxy environment of radio-loud quasars and radio-quiet quasars and concluded that there is 
no significant difference in the richness. Recent studies support this conclusion (e.g. Wold et 
al. 2001). The discrepancies between different studies may be caused partly by too small quasar 
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samples and by differences in sample selection of quasars. This situation will soon improve with 
the availability of a new generation of very large quasar surveys such as the 2dF quasar redshift 
survey (Groom et al. 2001) and the Sloan Digital Sky Survey (York et al. 2000). Although we do 
not deal with radio properties of quasars in this paper, our investigation of quasar environments 
will also provide a clue for understanding the radio character of quasar environments. 

The mean number of quasars per halo, {Nqso{M)) , and probability distribution of the 
number of galaxies around quasars, P(A''gai|A''Qso), used in this study can provide some use- 
ful features of the quasar environments. Furthermore, the spatial galaxy-quasar correlation 
function is used in order to quantify the galaxy environments around a quasar. Therefore, for 
the further investigation of environments and clustering of quasars and in order to constrain 
the quasar formation model, it is also necessary to predict spatial distribution of galaxies and 
quasars. We will show the results using the combination of cosmological A^-body simulation 
and SAM for formation of galaxy and quasar in the near future. 

We would like to thank T. T. Takeuchi for providing us with the reanalyzed data of 
the quasar luminosity functions derived from the 2dF 10k catalogue. We are also grateful to 
K. Okoshi, H. Yahagi and S. Yoshioka for useful comments and discussions. We also thank 
to the anonymous referee for a thorough reading of the manuscript and for his valuable sug- 
gestions and comments, which improved our paper very much. Numerical computations in 
this work were partly carried out at the Astronomical Data Analysis Center of the National 
Astronomical Observatory, Japan. This work has been supported in part by the Grant-in-Aid 
for the Scientific Research Funds (13640249) of the Ministry of Education, Culture, Sports, 
Science and Technology of Japan. 

Appendix. Star Formation and Gas Evolution 

In this appendix, we summarize our model of star formation and gas evolution. We use 
a simple instantaneous recycling approximation of model star formation, feedback and chemical 
enrichment. The following difference equations describe the evolution of the mass of cold gas 
Mcoid, hot gas Mhot, and long lived stars Mgtar at each time step. 

Mcoid = -M* + RM^ - (Al) 
Mhot = /5M„ (A2) 
M,t.r = M,-RM,, (A3) 

where M^, = Mcom/t* is star formation rate, R is the gas fraction returned by evolved stars, and 
(3 is the efficiency of reheating. In this paper, R = 0.25. The solutions of these equations are 
the following: 

t 



Meoid = M,^„idexp 



R + P) 



(A4) 



17 



Mhot = M°o, + /5AM„ (A5) 
Mstar = Ml, + (1 - R)AM,, (A6) 

where M^^i^,M^^^ and M^^^^. are the masses of cold gas, hot gas and long-hved stars from the 
previous time step, t is the time sine the start of the time step, and AM* = (M^^jj — Mcom) /(I — 
R + (3) is the mass of total formed stars. 

When a starburst occurs, stars are formed in a very short timescale. Thus, the starburst 
corresponds to r^/t — ^ in the above solutions. In this case, the changes of masses are given by 

Mcoid = 0, (A7) 

Mhot = M°„, + ^^^, (A8) 



M.ar = Ml, + il-§^ (A9) 



and the total star mass formed at starburst becomes 

AM,,burst = (AlO) 

i — iX + p 

From equation (AlO), we can obtain the mass of accreted cold gas onto a black hole (eq.[5]). 
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